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Abstract 

We argue that, at finite temperature, parity invariant non-compact elec- 
trodynamics with massive electrons in 2-1-1 dimensions can exist in both 
confined and deconfined phases. We show that an order parameter for the 
confinement-deconfinement phase transition is the Polyakov loop operator 
whose average measures the free energy of a test charge that is not an integral 
multiple of the electron charge. The effective field theory for the Polyakov 
loop operator is a 2-dimensional Euclidean scalar field theory with a global 
discrete symmetry Z, the additive group of the integers. We argue that the 
realization of this symmetry governs confinement and that the confinement- 
deconfinement phase transition is of Berezinskii-Kosterlitz-Thouless type. 
We compute the effective action to one-loop order and argue that when the 
electron mass m is much greater than the temperature T and dimensional 
coupling e^, the effective field theory is the Sine-Gordon model. In this limit, 
we estimate the critical temperature. Tent. = e^/87r(l — e^/127rm -|- . . .). 
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Gauge field theories in 2+1- dimensions have many interesting field theo- 
retical features such as gauge invariant local topological mass JI], H, fractional 
spin 1^ and statistics They have been of interest as model systems with 
somewhat more severe infrared divergences than their 3+1-dimensional rela- 
tives [0, 1^ and where one can study the dynamics of chiral symmetry break- 
ing 1^, 0, §]. Gauge field theories arise naturally in the description of lower di- 
mensional statistical models such as spin systems and the Hubbard model 
and an understanding of both their ground state and thermodynamic prop- 
erties are essential to the applications of mean field theory there. In this 
Letter, we shall show that parity invariant 2+1-dimensional quantum elec- 
trodynamics exhibits one more interesting property - when the electron has 
a large enough mass - it has a finite temperature confinement-deconfinement 
phase transition. 

The question of confinement has been extensively investigated in the more 
familiar context of quantized Yang-Mills theory at finite temperature and it 
is intimately related to the realization of a global symmetry involving the 



center of the gauge group |T0|, |TT[ . This symmetry transforms the Polyakov 
loop operator 

P{x) = trPexp ( drAoir, x)] (1) 



which measures gauge group holonomy in the periodic Matsubara time in 
the Euclidean path integral description of finite temperature gauge theory.^ 
Consider a gauge transformation 

v4/(r, x) = g'^{T, x)Af,{T, x)g{T, x) + ig'^ir, x)W ^g{T, x) (2) 

under the group SU{N) which has center Z^, the additive group of the in- 
tegers modulo N. In Eq.(|D g{T,x) can be periodic up to an element of the 
center of the group, g{l/T,x) = (yf(0, x)e^™/^. Under a gauge transforma- 
tion of the form 

P\x) = P(f)e2™/^ (3) 

Therefore, if the Z^r symmetry is not spontaneously broken, the correlators 
of Polyakov loop operators < P(xi) . . . P{xm)P\yi) ■ ■ ■ P\yn) >, are zero 

^We use units where Planck's constant, the speed of hght and Bohzniann's constant are 
one. For a discussion of the path integral formulation of finite-temperature gauge theory, 
see HI. 
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unless m = n modulo A^. The quantity 

F{xi,...,xm,yi,...,yn) = -rin(< P{xi)...P{x^m)P\yi)...P\yn) >) 

(4) 

is the free energy of the gluodynamic system in the presence of an external 
fundamental representation quark sources at positions xi,...Xm and anti- 
quark sources at yi, . . . ,yn- If the Zn symmetry is unbroken, the expectation 
value of a single loop, < P{x) >, vanishes and consequently the free energy 
F{x) of a single quark source is infinite. This is a signal of confinement - 
introducing a colored source into the confining system requires infinite energy. 
On the other hand, if the symmetry is spontaneously broken, F{x) is finite 
and characterizes the deconfined phase ||13|| . 

When dynamical quarks in the fundamental representation of the gauge 
group are present, the Polyakov loop operator does not characterize the con- 
fining phase, since, even if quarks are confined, screening can still take place. 
The coupling of gluons to quarks which are in the fundamental represen- 
tation, is invariant only under strictly periodic gauge transformations and 
therefore breaks the Zn symmetry explicitly. This is interpreted as the pos- 
sibility of fundamental quarks screening the color of an external source, so 
that the free energy of the source is always finite. 

In this Letter, we shall argue that the Polyakov loop operator can be used 
to study confinement in non-compact quantum electrodynamics even when 
dynamical electrons are present. In this case, the free energy of a distribution 
of external charges is 

with (anti-)periodic boundary conditions A^{l/T, x) = A^{0, x), ip{l/T, x) = 
— ?/'(0,x), ijj{l/T,x) = —ip{0,x). The gauge transformation 

A;{t,x) = A,{T,x)+V,x{r,x) , ij\r,x) = e*^^(^'^^V(r, x) , ij'{T,x) = ij{T,x)e-'^^'^^'^^ 

(6) 

is a symmetry of the action and measure if it preserves the (anti-)periodic 
boundary conditions, V /j^xi'^ / T , x) = V^x(0, ^) and x{'^/T,x) = x{^^^) + 
27m/ e. The coset of the group of all time- dependent gauge transformations 
modulo those which are periodic is the group Z, the additive group of the 
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integers. The Abelian Polyakov loop operator transforms under this global 
symmetry as 



exp yiYl'^i dTAo'{T, Xi)j = exp yiYl^^i dTAo{T, Xi)j exp ^ 

(7)^ 

Thus, if Z is not spontaneously broken, F{xi) defined by (|^) is infinite when 
the total charge of the external distribution is not an integral multiple of the 
electron charge, J^i^i 7^ integer ■ e. When the symmetry is broken, F{xi) 
can be finite. Thus, the nature of the realization of Z tests the ability of 
the electrodynamic system to screen charges which are not integral multi- 
ples of the electron charge and therefore probes the confining nature of the 
electromagnetic interaction. 

Z is the analog of the global Z^r symmetry of gluo-dynamics. However, 
in contrast to the case of a compact non-Abelian gauge theory, where Zjv 
is explicitly broken by dynamical quarks, and, due to compactness of the 
gauge group, charges which are not integer multiples of the quark charge are 
not available, the non-compactness of the gauge group of QED allows the Z 
symmetry to exist even in the presence of dynamical electrons. The electrons 
could be viewed as the analog of adjoint particles, either gluons or adjoint 
quarks, in QCD. 

At T = 0, and for the physical value of the electromagnetic coupling con- 
stant, 3+1- dimensional electrodynamics does not exhibit a confining phase. 
It is in the deconfined Coulomb phase at zero temperature and forms a De- 
bye plasma at finite temperature and density. There is a conjecture that, if 
the electron charge is increased so that e'^/Arc ~ 1, there is a phase transi- 
tion to a chiral symmetry breaking and perhaps confining phase [jn|. This 



phase, being in the strong coupling region, is difficult to analyze. In 1+1- 
dimensions, the massive Schwinger model is confining and the Z symmetry is 
not broken, at least when the temperature is much greater than the electron 
mass and the confinement scale is set by the dimensional electron charge e. 
On the other hand, it is known that when the electron mass is zero, the 



Z symmetry is spontaneously broken ||T5[. The symmetry breaking can be 
attributed to nonlocal effects of massless fermions. It can be argued that the 
phase transition between the broken Z and unbroken Z phases occurs for all 
temperatures at some value of the electron mass. 

One might also expect a confinement- deconfinement transition in the in- 
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termediate case of 2 + 1- dimensional electrodynamics. In that case, even at 
the classical level, the Coulomb potential is a marginally confining logarithm. 
Its entire spectrum is bound states, but the bound states can have arbitrarily 
large size. The free energy of a gas of charged particles is 

It is straightforward to compute the correlators of Polyakov loop operators 
(which are simply exponentials of the appropriate free energies, e~^'=' ''^). 
They have the scaling form 

i i<j 

with temperature dependent exponent, reminiscent of the spin-wave correla- 
tors in Gaussian spin wave theory in 2 dimensions |T^. 

It is interesting to ask how this result would be changed by radiative 
corrections and by thermal fiuctuations. This can be done by computing the 
effective action for the Polyakov loop operator. In electrodynamics at finite 
temperature, it is possible to use a gauge transformation to set the temporal 
component of the gauge field, Aq, independent of the Euclidean time. Then, 
the effective 2-dimensional field theory for Aq is obtained by integrating the 
other degrees of freedom from the path integral. What remains is an effective 



action for a static field Ao{x). In general, it was shown in ref. [jT^I that the 
effective action for the Polyakov loop operator in D+l-dimensional gauge 
theory is a D-dimensional sigma model with group-valued fields. In the 
present case of electrodynamics, the group-valued variables are 6*^"^°^^ and 
the effective field theory for Aq describes the appropriate sigma model. The 
Z symmetry is a periodicity of the effective action under the field translation 
Ao{x) ^ Aoix) + 27rT/e. 

The effective action obtained from integrating propagating fields from the 
path integral is non-local and non-polynomial in the remaining fields. It can 
only be regarded as a local field theory when the momenta of interest are 
much smaller than the masses of the fields which have been eliminated. In 
that case the effective action has a local expansion in powers of derivatives 
divided by masses. The effective action for Aq possesses such a local expan- 
sion. If the electron mass is sufficiently large that this expansion is accurate, 
the effective field theory for Aq can be approximated by a local field theory. 
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In 2+1-dimensions, the fermion mass operator constructed from the min- 
imal 2-component Dirac fermions is a pseudoscalar and therefore violates 
parity |TB|]. If included in the action, they can generate a parity violating 
topological mass for the photon by radiative corrections [0, 0. In this 
paper, we wish to study the case where the photon is massless. For this 
purpose, we study the model with Euclidean action 



S= d^x 



^F^^ + V^i(7 ■ (V + ieA) + m)^i + ^2(7 ■ (V + ieA) - m)V^2 



10) 



Where the parity transformation is the spacetime parity as well as ip'^ 2(2^') = 
7iV'2,i(a;), 4'i2{^') = — '^2,i(a;)7i with x' = {—xi,X2). The mass term ipiipi — 
ip2i^2 is a scalar. If m = in (|10D there is a 'chiral' symmetry under the 
transformation ipi{x) = Uijipj{x) where u G SU{2). The latter symmetry is 
broken at T = ||^, 0], however, since it is a continuous symmetry, it must be 
unbroken at any finite temperature in 2+1-dimensions. It has been argued 
that at finite temperature, the chiral transition is replaced by a Berezinskii- 
Kosterlitz-Thouless (BKT) transition [^]. In this paper, we shall consider 
the opposite limit of large mass, and show that there is indeed a BKT- 
transition corresponding to a confinement-deconfinement transition at some 
value of the temperature T. It is likely that this transition is in some way 
related to the chiral transition. 

At finite temperature, 2+1- dimensional QED contains three parameters 
with the dimension of mass, the electron mass m, the gauge coupling e^, and 
temperature T. The loop expansion is super-renormalizable |^, |^ and is an 
expansion in the dimensionless ratios e^/m and e^/T. We can compute the 
effective action for Aq(x) = a{x)\/T in a double expansion in the number 
of loops and in powers of derivatives of a{x). To order 1-loop and up to 
quadratic order in derivatives the effective action is 

S^s[Ao\ = I dx (^Z{m, ea/Vf)^Va ■ Va - V{m, ea/Vf)^ . (11) 

Here V is the effective potential for Aq arising from the fermion determinant 
and Z is obtained from expansion of the temporal components of the vacuum 
polarization function to linear order in — V^. To order one-loop, the effective 
potential is obtained from the fermion determinant in constant background 
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V{m,eAo/T) 



(Vol.; 



■ log det((-ii(9o - cAqY - + 



(12) 



where the fermions have anti-periodic boundary conditions in the direction. 

Dnsidering the ratio 

V2 + m2)/det(-92- 



The determinant can be computed by considering the ratio 

A(m, eAo/T) = det{{-ido - eAo)^ 
One finds 

sin2(eAo/2T) 



m 



(13) 



A(m,eAo/T)=n 



1 - 



cosh2(Afc/2r) 



= l[A^{m,eAo/T) 



(14) 



where A| = A;^ + are the positive eigenvalues of the operator — + m^. 
Eq.(|T^) holds in any dimensions. In 2+1-dimensions however one can perform 
the integral on k arising in log A(m, cAq/T), after taking the infinite volume 
limit. 



V{m,eAo/T) 
T2 



1 



(Vol. 



(2vr) 



logAg(m,eAo/T) 



TT 



^L^2(e-"^/^, eAo/T + tt) + Li3(e-"^/^, eAo/T + tt) 



T 



(15) 



where Li2{r,9) = — Jq dx\ia{l— 2x cos O+x"^) /2x and Li^i^r, 9) = Jq dxLi2{x,9)/x 
are the real parts of the dilogarithm and trilogarithm according to the con- 
vention of Ref. Eq. ([l5D shows the periodicity of the effective potential 
for eAo/T cAq/T + 27r. This is the residual gauge invariance. In 1 and 
3 dimensions the integral in k can only be performed for m = in which 
case it gives simple polynomial expressions. In the limit m = 0, the effective 



potential for Aq has been discussed in |22 



It is also straightforward to compute the term which contributes the lead- 
ing order in derivatives to the effective action, 



Z{m, ea/ VT) 



sinh m/T 



m 



-m/T 



VlTim [^coshm/T + cosea/Vr ^ (^cosh m/T + cos ea/^T 

(16) 
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The critical behavior of the 2- dimensional model defined by Eqs.(|Tl|), (p^) 



and (IC), can be understood by comparing it with the sine-Gordon model 
in two- dimensions. That this comparison can be reliably performed can be 
seen by the study of the harmonic content of (pISl). 



^ / r}m\ 

\/(m,ea/VT) = -— J^e""'"/^ 1 + ^ cos(n(ea/Vr + 7r)) . (17) 



2 oo 



^ ntl V T J 

Consider then the large m limit, T/m and e^/m small with finite e^/T. In 
this limit, the higher harmonics are small perturbations to the potential 

y(m,ea/yT) = — e-'"/^cos(ea/Vr) , (18) 

71 

which is the sine-Gordon potential. Amit et al. showed that in the 

sine-Gordon model any perturbation of the type cos(n/30) to a sine-Gordon 
potential a cos(/50) / are irrelevant for the critical behavior of the model. 
By analogy with the spin wave plus Coulomb gas model, it was also proven 
in Refs. that a critical line for a BKT |2^, ^ phase transition in the 
sine-Gordon model with a logarithmic potential starts at the point (a,/3^) = 
(0,87r). We can then conclude that also in 2-l-l-QED at finite temperature 
there is a BKT phase transition, with a critical line in the [m/T, e^/T) plane 
starting at {m/T,e^ /T) = (cxD,87r). The critical temperature for this transi- 
tion (up to 1-loop order) can be computed from Eq.(0) and Eq.([T^) as 




e2 



127rm 



+ ... . (19) 



This is the critical value of the coupling constant originally found by Coleman 
in his discussion of bosonization of the massive Thirring model . 

Note that the vacuum expectation value of Aq in the deconfined phase, 
where the Z symmetry is spontaneously broken, is 

27mT , . 

<Ao>=^— . (20) 

In a semiclassical analysis, this expectation value contributes an imaginary 
chemical potential for the electron action. However, this chemical potential 
can be absorbed by shifting the Matsubara frequency. Thus, the semiclassical 
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thermodynamics do not suffer from the difficuhies of the meta-stable Zn 
phases of QCD l^, |3ll . 
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